Near-field EM wave scattering from random self-affine fractal metal surfaces: spectral dependence 
of local field enhancements and their statistics in connection with SERS 
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By means of rigorous numerical simulation calculations based on the Green's theorem integral equation formu- 
lation, we study the near EM field in the vicinity of very rough, one-dimensional self-affine fractal surfaces of 
Ag, Au, and Cu (for both vacuum and water propagating media) illuminated by a p polarized field. Strongly lo- 
calized enhanced optical excitations (hot spots) are found, with electric field intensity enhancements of close to 
4 orders of magnitude the incident one, and widths below a tenth of the incoming wavelength. These effects are 
produced by roughness-induced surface-plasmon polariton excitation. We study the characteristics of these op- 
tical excitations as well as other properties of the surface electromagnetic field, such as its statistics (probability 
density function, average and fluctuations), and their dependence on the excitation spectrum (in the visible and 
near infrared). Our study is relevant to the use of such self-affine fractals as surface-enhanced Raman scattering 
substrates, where large local and average field enhancements are desired. 
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I. INTRODUCTION 

Recently, single molecule probing by means of Surface En- 
hanced Raman Scattering (SERS) has been reported both on 
Ag single nanoparticleslil and on Ag colloidal aggregates^ the 
latter exploits the extremely large near infrared (NIR) Raman 
scattering cross sections of dye molecules.0 Bearing in mind 
how inefficient normal spontaneous Raman scattering is, en- 
hancement factors of 10 14 , or larger are required to achieve 
single molecule detection. El 

Typically, SERS is known to yield Raman signals enhanced 
by a factor Q ^ n lfl 6 with respect to those of conventional Ra- 
man scattering.a~Q Two mechanisms are responsible for such 
enhancement factors: the surface-roughness-induced intensi- 
fication of the electromagnetic (EM) field both at the pump 
frequency and at the Raman-shifted frequency (EM mecha- 
nism), and the charge-transfer mechanism. The former mech- 
anism is widely accepted to be the most relevant one from 
the quantitative standpoint, providing gains of Gem > 10 4 
in most experimental configurations. Extensive theoretical 
work has been devoted to the explanation of the EM mech- 
anism (cf. e.g. the reviews Refs. and the consensus 
is that what underlies such EM field enhancement (FE) fac- 
tors is the-,roughness-induced excitation of surface-plasmon 
polaritonscl (SPP), either propagating along a continuous sur- 
face (extended SPP), confined within metal particles (particle 
plasmon resonances), or even confined due to Anderson local- 
ization (localized SPP). 

In light of the SERS enhancement factors estimated for sin- 
gle molecule detection it is evident that, in addition to the well 
known average SERS enhancement factors, extremely large 
EM fields must appear in the vicinity of SERS substrates, even 
if the charge-transfer mechanism is also known to be espe- 



cially intense, as in Ref. |[ This is supported by the observa- 
tion, through photon scanning tunneling microscopy (PSTM), 
of very intense and narrow EM modes (called hot spots) on. 
rough metal surfacesE3~E3 or rough metal-dielectric films.E3 
Interestingly, these rough metal surfaces used as SERS sub- 
strates posses, in some cases, scaling properties within a suf- 
ficiently wide range of scales (physical fractality). The sub- 
strates can presenL&elf-similarity, as the widely employed col- 
loidal aggregatesJl3'Ejil3 or self-affinity, as in the case of de- 
posited collnids-or evaporated or etched rough surfaces and 
thin filmse-ErEl 

Therefore, inasmuch as the quantitative evaluation of the 
surface EM field is central to the SERS effect, knowl- 
edge of the EM scattering process for surface models as 
realistic as possible is obviously needed. In recent years, 
the theoretical efforts have been directed towards either 
describine-jfliflMjig-h! approximate methods realistic surface 



model sB 




the full EM theory to study sim- 
though introducing increasingly 



plistic surface mode 
complex properties.E3 

In this paper we study the near EM field scattered in the 
vicinity of rough, one-dimensional self-affine fractal surfaces 
of Ag, Au, and Cu, with the aim of determining the appear- 
ance of strong local optical excitations (hot spots) and charac- 
terizing them with regard to their spatial and spectral width, 
their polarization, and their excitation spectra; in addition to 
that, the global optical response of such fractal surfaces will 
be studied through the statistical properties of the surface EM 
fields. Both local and global responses are discussed in light 
of the influence on SERS. For this purpose, we make use of 
numerical simulation calculations basecLon the Green's the- 
orem integral equation formulation£frEj rigorous from the 
classical EM standpoint. Unlike recent work also for self- 
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affine fractals (though sensitively/Jess rougher) based strictly 
on magnetic field calculations,! 2 ^!^ the magnitude that natu- 
rally arises in this formulation when applied to ID surfaces 
and p polarization (the one relevant for its light-SPP coupling 
selectivity), we fully characterize here the surface and near 
electric field components (crucial in SERS and other nonlin- 
ear optical effects) through simple expressions in terms of the 
magnetic field and its normal derivative. The necessary details 
of the theoretical formulation are given in Sec. £§ The local 
optical excitations are studied in Sec. |l| and their statistical 
properties in Sec. |v| leaving for Sec.^V] the conclusions of 
this work. 



II. SCATTERING FORMULATION 
A. Surface integral equations 

The scattering geometry is depicted in Fig. [[]. A rough 
metal surface z — ((x) is the substrate onto which 
molecules are adsorbed in SERS typical experimental con- 
figurations. The semi-infinite metal occupying the lower 
half-space [z < ((x)] is characterized by an isotropic, ho- 
mogeneous, frequency-dependent dielectric function e < (u>). 
From the medium of incidence, characterized by a frequency- 
dependent dielectric function e > (oj), a monochromatic, lin- 
early polarized incident beam of frequency us impinges on the 
interface at an angle 0o, measured counterclockwise with re- 
spect to the positive z axis. The polarization is defined as 
shown in Fig. |IJ: the magnetic (respectively, electric) field 
is perpendicular to the xz plane for p (respectively, s) polar- 
ization, also known as the transverse magnetic (respectively, 
transverse electric) one. 

We restrict the analysis to ID surfaces (invariant along the y 
direction), without any loss of generality as far as the.physics 
underlying the SERS EM mechanism is concerned.EjE3 It has 
been shownETKlJ that the latter condition simplifies consider- 
ably the formulation based on the integral equations resulting 
from the application of Green's second integral theorem (with 
the help of the Sommerfeld radiation condition). In such cir- 
cumstances, the starting 3D vectorial problem can be cast into 
a 2D scalar one, where the unknown is the y-component of 
either the magnetic field [Hy (r, uj), with r = (x, z)] for p 
polarization, or the electric field [E y s \r,uj)] for s polariza- 
tion; no depolarization takes place when pure s or p polarized 
fields are incident. This simplification is very convenient from 
the analytical and numerical point of view, r aii< | LYields straight- 
forwardly the far field scattered intensity .EjlEI We point out 
that, if the electric (respectively, magnetic) field is needed for 
p (respectively, s) polarization, it can be calculated on the ba- 
sis of Maxwell equations, as we shall see below. 

Let us focus on the electric field calculation for the case 
of p polarization. Evidently this is the most relevant one for 
our SERS problem (only p polarized light can excite in the 
present configuration the SPP responsible for the EM field en- 
hancements), and also to other interesting problems such as 
second harmonic generation on metal surfacesaEH or PSTM 



studiesOTl 2 ] As mentioned above, the integral equation for- 
mulation is simplified if written in terms of the magnetic field 
amplitude. Our monochromatic incident field ofJiequency uj 
is a Gaussian beam of half-width W in the formra 
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where k e = n^uj/c and n£ = Ve>. From now on, since a 
time harmonic dependence e~ lut is assumed, the functional 
dependence on frequency will be omitted unless necessary for 
the sake of clarity. The surface integral equations that fully 
describe the EM linear scattering problem for p polarization, 
in the geometry of Fig. [j], are 
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where H, 



; ' y \r) and Hy P,< \r) are the magnetic field in the 



upper (z > Q and lower (z < Q semi-infinite half-spaces, 
and the normal derivative is defined as d/dn = (n • V), with 
n = 7- 1 (-C'(z),0,l) and7 = (l + (C'(x)) 2 ) 1 / 2 . The 2D 
Green's function G is given by the zeroth-order Hankel func- 



tion of the first kind H, 
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The four integral equations (2.2) fully describe the scatter- 



ing problem for p polarization in terms of the y component of 
the magnetic field. Analogous integral equations can be ob- 
tained for s polarization dealing with the y component of the 
electric field. In order to solve for the surface field and its nor- 
mal derivative, defined as the functions H(x) and j~ 1 L(x), 
two of the integral eq uatio ns (no te tha t they are not indepen- 
dent), typically Eqs. ( 2.2a ) and ( 2.2c ), are used as extended 
boundary conditions, leading to two coupled integral equa- 
tions once one invokes the continuity conditions across the 
interface: 

H(x) = ff(»(r) |, =c c+) ((e) = H<\v) U c c- )(x) , (2.3a) 



j^ 1 L(x) 



dn 
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dn 
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with ^^'(x) = lim £ _>o (C( x ) i £ )- Th e resulting system of 
integral equations can be numerically solved upon convert- 
ing it into aj-Sjystem of linear equations through a quadra- 
ture scheme,E!l t he un knowns being H(x) and L(x). Then 
Eqs. (2.2a) and (2.2d) permit to calculate the scattered mag- 
netic field in the upper incident medium and inside the metal, 
respectively. 

But what if the magnitude of interest is the electric field? 
This is indeed the situation in SERS where the surface elec- 
tric field locally excites the molecule vibrations that produce 
the Raman-shifted radiation that is detected. In Refs. p^ , p7[ , 
the EM field enhancement factor has been defined as the nor- 
malized magnetic field intensity: 



a H (uj) 



H, 
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(2.4) 



Even if the enhancement factor thus defined closely resem- 
bles the correct total electric field enhancement factor, we are 
evidently losing information about the different electric field 
components, in turn relevant to the SERS polarization selec- 
tivity. 



B. p polarization: Electric field 

In order to obtain the electric field components from the 
y component of the magnetic field, use can be made of the 
Maxwell equation 
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In the incident medium, Eq. (2.2a) provides the only non-zero 
com ponent of the magnetic field. Use of Maxwell's equation 
(2.5) leads to the following electric field components: 
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These equations can be rewritten in terms of the source func- 
tions H(x) and L(x) as follows: 
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where the explicit form of the Green's function has been taken 
into account, leading to the appearance of 1st and 2nd order 

Hankel functions of the first ki nd h \ x \ H% ■ F° r tne Gaus- 
sian incident field given by Eq. (2.1 ), the electric field compo- 
nents are 

M^)(r) = -i-fffe^r) \ icos8 {l + w{x,z)) 
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(2.8b) 
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Equations (2.7) and (2.8) provide the electric field compo- 
nents in the incident medium of the resulting p-polarized EM 
field, incident plus scattered from the rough surface. The scat- 
tered electric field involves an additional surface integral in 
terms of the source functions, previously obtained (numeri- 
cally) from the above mentioned coupled integral equations. 
Analogous expressions, not shown here, for the correspond- 
ing electri c field components inside the metal can be obtained 
from Eq. ( 2. 2d ). On the other hand, recall that a similar pro- 
cedure can be straightforwardly developed to yield the mag- 
netic field components in the case of s-polarized EM waves as 
surface integrals in terms of the surface electric field (y com- 
ponent) and its normal derivative. 



C. p polarization: Normal and tangential surface electric field 

It should be pointed out that when trying to evaluate the 
elect ric field on the surface, or even very close to it, from Eqs. 
(2.7), non-integrable singularities appear associated with the 
Green's functions derivatives for vanishing arguments. Use of 



3 



expressions ( 2.7) for the evaluation of the electric field close to 
the surface will produce unphysical results. To deal properly 
with this situation, more care s hould have been taken in doing 
the derivatives of the integral ( 2.2a ) describing the magnetic 
field, whose, integrand already exhibits singularities, though 
integrable] 2 JrE] A simple way to work around this problem 
consists of evaluating the electric field at the surface itself, and 
we have found very simple relations connecting the normal 
and tangential components of the electric field (see Fig. [j]) 
with the surface magnetic field and its normal derivative [cf. 
Eqs. (2.3)]. These are 



ic _ 1 dH(x) 
7 



dx 



E 



(x) = 



ic 



-— 7 - i i(x). 



(2.9a) 
(2.9b) 



These expressions are extremely useful, for they facilitate 
considerably our study of the SERS EM mechanism. 

Thus, taking advantage of the one-dimensional scattering 
geometry (which, although is not general., is not inappropri- 
ate to study the SERS EM mechanismc§E3), we only have 
to deal with the y-component of the magnetic field to obtain 
a simplified, basically exact solution to the scattering prob- 
lem from the classical EM viewpoint, at the excitation fre- 
quency. The numerical solution of the resulting integral equa- 
tions yields the surface magnetic field and its normal deriva- 
tive as the main results. The drawback of working with the 
magnetic field, when the quanti ty o f interest is the electric 
field, is avoided by expressions (23), which allows us to ob- 
tain the surface electric field with the only additional algebra 
of calculating a spatial derivative. 

We now properly define the electric field enhancement fac- 
tors for either a single component or the total field as the nor- 
malized intensities: 
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D. Numerical implementation 



The numerical procedure has been implicitly outlined 
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Self-affine 



above; further details have been given in Ref. 
random fractal surfaces numerically generalad-by means of 
Voss' fractional Brownian motion algorithmoS are studied. 
This kind of fractals exhibit self-affine scaling properties in 
a broad spatial range£3 and have properties that resemble 
those of some SERS substrates—such as cold-deposited metal 
films or etched metal surfaces.E]i£j In the numerical calcula- 
tions, surface realizations of length L = 10.29/im, consist- 
ing of N p = ri,* iV sampling points obtained by introducing 
n, = 4,6,8 or 10 cubic-splined interpolating points into a 
sequence of N — 201 points extracted from each generated 



fractal profile with Nf — 1024 points; note the considerably 
larger sample density with respect to that of Ref. ^7|. The sta- 
tistical properties of the physical quantities of interest will be 
calculated on the basis of Monte Carlo simulations for an en- 
semble of fractal realizations. 



III. LOCAL FIELD ENHANCEMENT: HOT SPOTS 

We now turn to the investigation of the occurrence of very 
large near EM field enhancements. Particularly, we will con- 
centrate on self-affine fractals with Hurst exponent H = 0.1 
(namely, local fractal dimension Dj = 2 — H = 1.9), which 
have been shown in Ref. ^7] to give rise to large surface mag- 
netic fields. The lower scale cutoff £x, ~ Wif-H 11 has been 
chosen to resemble that of SERS substtatesJiSEI The upper 
scale cutoff, typically ~ 50/^m, is considerably larger than 
the illuminated area L, and this is in turn sufficiently (in order 
to avoid finite length effects) larger than the incoming wave- 
length (0.4 ^m< A < 1.3fim). Thus physical scaling is mean- 
ingful for the relevant interval of this scattering problem. The 
effect of further reducing the lower scale cutoff will be inves- 
tigated elsewhere;t3 in this regard, it should be recalled that 
the minimum scale relevant to the far field pattern has been 
studied for Koch fractals.Ea 



A. Near Field Intensity Maps 

In Fig. |[ the intensity (on a logarithmic scale) of the elec- 
tric and magnetic near field in the vicinity of a self-affine Ag 
surface with D = 1.9 and rms height 5 = 514.5 nm, in a 
particular region (of about 1 x 1/im 2 ), is shown for normal 
incidence with light of wavelength A = 514.5/im; in addition, 
the intensities of the two different components (x and z) of 
the electric field are separately shown. Before analyzing the 
results, some comments are in order with regard to the nu- 
merical calculations. Whereas the electric and magnetic fields 
in vacuum away fro m the interface are given by the integral 
equations (2.7) and (2.2a), respectively, and similarly for the 
EM field inside silver, their corresponding values on the inter- 
face are d irectly obtained from the source functions through 
Eqs. (O) and ( p|) As mentioned in the preceding section, 
Eqs. (2.2a) and (2.7) exhibit singularities upon approaching 



the surface, so that they are not accurate at points very close 
to the surface, typically within distances smaller than the sur- 
face sampling interval. Thus the EM field intensity in Fig. || 
at distances from the surface smaller than AL/N p for the mag- 
netic field and 6L/N p for the electric field are obtained from 
the weighted values on the two closest sampling points on the 
surface, taking explicitly into account for poin ts in side silver 
the continuity conditions for the magnetic [Eq. (2.3)] and elec- 
tric field components. Despite that, some slight (not inherent) 
mismatch might still appear when entering into the surface 
field area, mostly in the intensities of the electric field compo- 
nents. 
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It is, of course, expected that such continuity conditions 
across the interface could be roughly observed in the calcu- 
lations even if we were not considering the above explicit 
matching. On the one hand, the continuity of the tangential 
component of the magnetic field is neatly seen in Fig. Eld); on 
the other hand, the continuity of the tangential electric field is 
appreciable in Fig. ^(b) and (c) through the continuity of the 
x (respectively, z) component of the electric field at locally 
flat (respectively, vertical) parts of the rough surface, whereas 
the discontinuity of the normal component of the electric field 
(continuity of the normal component of the displacement vec- 
tor) is inferred from the discontinuity of the x (respectively, 
z) component of the electric field at locally vertical (respec- 
tively, flat) areas. Incidentally, note also that the EM field 
inside silver decays very rapidly as expected from the Ag skin 
depth d = (c/cj)(-e < )- 1 / 2 « 27 nm (cf. Ref. || for the Ag 
dielectric constant). 

It is evident from Fig. ^ that the maximum local EM fields 
are located right on top of the Ag surface, whereupon some 
particularly bright spots appear. Thus we next plot in Fig. || 
the intensities of the surface EM fields (including electric tan- 
gential and normal components) for the surface area shown in 
Fig. ||, including the surface profile. Very narrow peaks sur- 
rounded by dark areas are observed, whose widths are well be- 
low half wavelength of the SPP (X S pp/2 = 243 nm). Some 
of these peaks can be considered as optical excitations (hot 
spots) where very large local FE's occur, such that nonlin- 
ear optical processes would be strongly enhanced thereinE In 
particular, the largest in Fig. |3|, shown in the inset, is of the 
order of a « 6 • 10 2 . Note that our calculations identify the 
local electric field component that is responsible for such FE: 
the normal component. In Fig. |], the near electric field in the 
vicinity of this hot spot (zooming in Fig. ^(a)) is given; inter- 
estingly, it is associated with a surface peak. 

It has been experimentally shown by near-field microscopy 
that such optical excitations rapidly disappear upon changing 
the frequency of the incident radiation.E3 Our rigorous calcu- 
lations corroborate those experimental observations, as seen 
in Fig. ||, where the surface electric field intensity is plotted 
for several incident wavelengths close to A = 514.5 nm. For 
the wavelengths A =495.9 and 539.1 nm, the high intensity 
spot is still clearly visible, though less bright. It fades away, 
however, for larger frequency shifts, becoming barely visible 
for A > 563.6 nm or A < 476.9 nm (about 10% frequency 
shift). 

The hot spot shown in Fig. ^| is strongly polarized along 
the normal to the surface. This is extremely relevant to SERS 
spectroscopy, since it might impose selection rules to the vi- 
brational modes of the adsorbed molecule. Is it possible to 
find hot spots with different polarizations? Only in the case of 
silver at wavelengths close to the surface plasma wavelength, 
we have found certain spots strongly polarized along the tan- 
gential direction too, though weaker than those normally po- 
larized. These tangentially polarized hot spots exhibit FE fac- 
tors not larger than a t ~ 10 2 , and are typically located in 
regions presenting larger values of a n : this will be discussed 
elsewhererJ 

The occurrence of local optical excitations has been studied 



for the same self-affine surface profile illuminated with differ- 
ent excitation wavelengths, and also for metals such as Au and 
Cu. Although not shown here, similar normally polarized hot 
spots are found on a broad spectral range on all the self-affine 
surfaces of Ag, Au, and Cu. We now discuss some of the 
characteristics of thespiiot spots. 

It has been arguedEil that these hot spots are due to An- 
derson localization of SPP. Theoretical works based on a 
dipolar model, also demonstrate the possibility of creating 
strongly confined and intense excitations on self-affine frac- 
tal surfaces;E3 in fact, in the case of random metal-dielectric 
films, Arijderson localization of surface plasmon modes is 
predicted.EJ Our numerical calculations, not subject to dipo- 
lar (and quasi-static) restrictions, indeed reveal the existence 
of these kind of optical excitations and, although compatible 
with the possibility of them being due to Anderson localiza- 
tion of SPP, do not permit to draw further conclusions in this 
respect. Our scattering geometry involving the interaction be- 
tween a propagating beam of light and a metal surface does 
not lend itself well for the characterization of the SPP Ander- 
son localization phenomenon. To that end, the study of the 
propagation and transmission of SPP through rough surfaces 
would be more adequate.EZI Only the fact that light can couple 
into these, possibly localized, SPP modes through the rough- 
ness can be inferred from our calculations and from the typical 
PSTM configurations^ 

It is also worth pointing out that the extremely rapid decay 
and widening of the optical excitations in the near field (see 
Figs. U and implies that PSTM images taken at certain dis- 
tance from the surface, leaving aside the rounding effects of 
the tip, will manifest themselves as much wider pacLweaker 
optical excitations, and this seems to be the case.ErE 2 ] Direct 
probing of hot spots, on the other hand, could be carried ouLhji 
the nonlinear effects of physi- or chemi-sorbed molecules.B tfl 

B. Spectral dependence 

In order to analyze the polarization and spectral dependen- 
cies of the optical excitations, we present in Fig. ^ the max- 
imum local FE values found at Ag, Au, and Cu self-affine 
surfaces with D = 1.9 and 5 — 514.5 nm, obtained from nu- 
merical calculations of the surface EM field for an ensemble 
of N r = 60 realizations generated as mentioned above (only 
the data from the central half of each realization are used). 
The results for weaker self-affine surfaces (S = 102.9 nm) 
used in Ref. ^7] to compute magnetic FE's are also shown. In 
addition, the case of having water as incident medium (sol- 
vent) has been analyzed, though in Fig. ^ only the results for 
H20/Ag are shown. Several remarks are in order with regard 
to Fig. | 

Very large FE's appear for a wide spectral range covering 
the visible and entering into the NIR. In the red and NIR parts 
of the spectrum, all three metals being studied behave sim- 
ilarly, giving rise to hot spots exhibiting strongly enhanced 
electric field intensities normal to the surface, the tangential 
component tending to vanish. This behavior can be under- 
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stood in accordance with the-spectral evolution of the Ag, Au, 
and Cu dielectric constants,!^ all showing increasingly large 
negative real parts, tending to the perfectly conducting limit 
e — > — oo that predicts vanishing tangential electric fields. For 
wavelengths A ^> 1240 nm, however, FE's are expected to 
slowly decrease as the surface is "seen" by the incoming ra- 
diation of increasing wavelength as increasingly flatter. Other 
calculations, not shown here, indicate that this is the case for 
A > 2/xm. In fact, this decay is observed at lower wavelengths 
(within the spectral interval covered by Fig. ||) for the fractal 
surface with 5 = 102.9 nm. 

The optical responses of Au and Cu manifest significant dif- 
ferences with respect to that of Ag in the blue part of the spec- 
trum. The onset of interband transitions, which takes place in 
Ag at A » 300 nm unlike in Au and Cu (slightly below 600 
nm), makes the difference inasmuch as such transitions con- 
stitute a strong absorption mechanism. Consequently, FE's 
should be significantly reduced for wavelengths below the on- 
set threshold, as is evident in Fig. |^ for Au and Cu below 
A w 600 nm, but not seen for Ag since the lower wavelength 
considered in Fig. ^ is above the Ag onset threshold. More- 
over, silver surfaces at small incoming wavelengths approach- 
ing the surface plasmon wavelength (but above the onset of 
interband transitions) present strong local optical excitations 
tangentially polarized, as mentioned above. These tangential- 
electric hot spots can lead to local FE nearly comparable to 
those corresponding to the normal-electric ones at such wave- 
lengths, although more than an order of magnitude weaker 
than those obtained at larger wavelengths [see Fig. ^|(b) and 
(c)]. This has important implications in SERS, since Ag sub- 
strates of the kind studied here, when illuminated at wave- 
lengths A < 600 nm, could enhance the Raman signal coming 
from a vibrational mode of the molecule sensitive to the tan- 
gential electric field, as well as those sensitive to the normal 
electric field (typically established as predominant according 
to SERS selection rulesQ). 

Finally, note that using water as solvent does not introduce 
significant changes in the qualitative and quantitative behav- 
ior of the maximum local FE. In addition, we would like to 
point out that the magnetic FE follows qualitatively (and al- 
most quantitatively) the normal electric FE. 



IV. STATISTICAL PROPERTIES OF THE SURFACE FE 

In this section we study the statistical properties of the FE's 
occurring on the surface of self-affine fractal profiles with 
fractal dimension and rms roughness as in the preceding sec- 
tion. These properties are obtained from Monte Carlo numer- 



ical simulations results performed as described in Sec. IIP 



Preliminary results based on magnetic field calculations for 
weaker fractal surfaces have been presented in Ref. E7l 



A. Probability Density Function 

In Sec. |ni| we have found, by direct observation of the cal- 
culated near-field excited in the neighborhood of the interface, 
that very large fields can be excited at the surface. These en- 
hanced excitations are coupled to the surface by the surface 
roughness and, thus, depend strongly on its properties. The 
question then arises as to how probable these values are and, 
in turn, how the EM field intensity is distributed over the frac- 
tal surface. In Fig. ^, we show the PDF of the surface EM 
FE (including separately tangential and normal components) 
for self-affine fractal surfaces with D — 1.9 and 6 =514.5 
nm for two incoming wavelengths A = 514.5 and 1064 nm, 
averaging over various angles of incidence. For the sake of 
comparison, the results for fractal surfaces with smaller rms 
height S =102.9 nm, and also with both smaller fractal di- 
mension D = 1.2 and rms height 6 =102.9 nm are included 
(the latter as used in Ref. for magnetic FE calculations). 

For the smoother surface, the resulting PDF is a narrow 
distribution centered at the surface EM FE value for a flat 
metal surface, as expected for such a weakly rough surface 
and in agreement with Refs. wherein similar results 

were shown for the intensity of the magnetic field (typically, 
0ff ~ \1 + R\ 2 , R being the corresponding Fresnel reflec- 
tion coefficient). Recall that the electric field components on 
flat metal surfaces follow a t ~ cos 2 0o|l + R\ 2 and a n ~ 
sin 2 1 1 - R\ 2 , so that the PDF distributions for D = 1.2 in 
Fig. [7] are only significant for FE's approximately in between 
the minimum and maximum expected values for the different 

For rougher surfaces, however, the surface EM no longer 
resembles the flat surface result, presenting alternating dark 
and bright regions, and giving rise with increasing roughness 
parameters to very bright hot spots surrounded by large dark 
regions. The corresponding PDF becomes wider, turning into 
a slowly decaying function that is maximum at zero and ex- 
hibits a long tail for large FE values (see Ref. |27| for the mag- 
netic FE PDF for the self-affine fractal with D — 1.9 and 
6 =102.9 nm). Upon comparing Figs. [7] (a) and (b), it is ev- 
ident that moderately large er t 's become feasible at A =514.5 
nm as well as very large <x n 's, whereas only the er„'s are ex- 
pected to be intense at A =1064 nm. 



B. Average and Fluctuations 

As a result of the change in the surface EM field PDF for 
increasing surface roughness parameters, the moments of the 
distribution are also modified. Particularly relevant are the av- 
erage and the statistics of the fluctuations, as they can also 
provide some information about the global response of larger 
surfaces (of the order of centimeters) under broad beam illu- 
mination. 

In Fig. [| we present the spectral dependence of the mean 
FE for the same self-affine fractal surfaces whose local FE 
where shown in Fig. ^ In fact, the qualitative behavior of 
the mean FE does not differ substantially from that exhibited 
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by the maximum local FE. Basically, there is a broad excita- 
tion spectra for the rougher fractals (slightly narrower for the 
smoother fractal), covering the visible and the NIR (at least up 
to A = 2/im). The behavior is similar for Ag, Au, and Cu and, 
in all cases, the electric field is predominantly normal to the 
surface. The blue part of the excitation spectra reveals, on the 
other hand, a rapid decrease for Au and Cu associated with the 
onset of interband transitions, whereas large normal FE's can 
still be found in that spectral region for Ag fractals as well as 
an increase of the tangential electric field upon approaching 
the surface plasma wavelength (but still above the threshold 
wavelength of interband transitions). The latter blue tangen- 
tial electric FE increase is slightly larger when water rather 
than air constitutes the propagating medium. On the other 
hand, it should be emphasized that the rougher fractal sur- 
faces used in Fig. || give rise to an estimated SERS FE factor 
(Gem ~ 10 5 , in good agreement with the phenomenological 
factor experimentally induced.B 

The absorption spectra are shown in Fig. || for the sake of 
comparison. It is evident that the absorption spectrum does 
not resemble the qualitative behavior of the excitation spec- 
trum in Fig. m Therefore, for this kind of fractal surfaces 
yielding wide excitation and absorption spectra, the maximum 
absorption region as experimentally obtained from the absorp- 
tion spectra cannot be straightforwardly related to the opti- 
mum excitation wavelength. Furthermore, it should be em- 
phasized that strong absorption can even be associated with 
very low surface EM fields (and thus the substrates being 
SERS inactive), as is the case of Au and Cu self-affine fractals 
in the blue spectral region. In other substrate configurations, 
nonetheless, the contrary might be the case, and absorption 
bands can be used to identify excitation bands resulting in 
strong surface FE (substrates becoming SERS active), as in 
rough surfaces presenting surface shape plasmon resonances, 
such as colloidalraggregates,Erla or in gratings diffracting into 
propagating SPPB In summary, one has to carefully interpret 
absorption spectra when using such information to determine 
the appropriate SERS (or any other surface optical nonlinear 
effect) excitation frequency. 

Finally, we show in Fig. [l(] the spectral dependence of the 
FE fluctuations. In accordance with the previously discussed 
FE PDF widening with increasing surface roughness, it is ob- 
vious that the rougher the surface, the larger the fluctuations. 
And for sufficiently rough surfaces, the fluctuations can be 
even larger than the average, as seen upon comparing Fig. 
with Fig. p| Actually, the FE fluctuations rather than the aver- 
age provide a good estimate of how probable and how bright 
hot spots are. Indeed, the spectral dependence of the FE fluc- 
tuations in Fig. [l(] closely follows that of the maximum local 
FE shown above in Fig. 



V. CONCLUSIONS 

By means of a rigorous Green's theorem integral equation 
formalism, we have studied the occurrence of strong local op- 
tical excitations (hot spots) on self-affine fractal surfaces of 



Ag, Au, and Cu. The statistics of the surface field fluctuations 
that produce these strong excitations have also been studied. 
The formalism exploits the scalar character of the resulting 
integral equations for one-dimensional surfaces illuminated 
with linearly s- or p-polarized light, by treating the problem 
in terms of the electric or magnetic field, respectively. In the 
case of p polarization, which is the relevant one in our prob- 
lem due to the SPP excitation selectivity, we have calculated 
the electric field from the resulting only nonzero component 
of the magnetic field and its normal derivative on the surface. 
The problem is studied numerically by means of Monte Carlo 
simulations of the interaction of light with self-affine metal 
fractals whose profiles were obtained from the trace of a frac- 
tional Brownian motion. The appearance of hot spots and their 
statistics have been determined for a broad spectral range of 
the incoming light (400 nm < A < 1300 nm). 

We have found hot spots on self-affine fractals with frac- 
tal dimension D = 1.9 and rms height <5 =514.5 nm. These 
hot spots constitute very strong and narrow (considerably nar- 
rower than half of the SPP wavelength) surface EM field ex- 
citations, with very selective excitation spectra (both tem- 
porally and spatially). Typically, they give rise to local FE 
strongly polarized along the normal to the rough surface. The 
largest ones we have found yield local SERS FE factors of 
Gem ~ 10 7 , and appear for a wide range of incoming wave- 
lengths covering the visible and NIR up to A w 2fim. In- 
terestingly, our results reveal that weaker, tangentially polar- 
ized hot spots can be found in Ag fractals for small excitation 
wavelengths (blue or smaller). 

All these features are not incompatible with the suggestion 
that Anderson localization of SPP is the underlying physical 
mechanism responsible for such optical excitations (e.g. the 
exponential decay of the SPP transmission versus rough sur- 
face lengthpCpuld certainly provide some indication of such 
mechanismO), although no direct evidence of this can be ob- 
tained from our scattering geometry and calculations. 

The PDF of the surface EM field for those self-affine frac- 
tals exhibiting hot spots is a slowly decaying function with a 
significant tail for large surface FE's. It differs substantially 
from that for smooth surfaces for which the PDF is a narrow 
distribution centered at the value of the EM field expected on 
a flat metal surface. 

The mean FE acquires considerable values in a broad spec- 
tral region. For the three metals considered, in most of the 
visible and NIR (A < 2/j.m) excitation regions, the component 
responsible for this enhancement is the electric field normal 
to the surface. We have found that, for these self-affine fractal 
surfaces the average SERS FE factors are (Gem) ~ 10 5 . 

We have also analyzed the spectral dependence of the sur- 
face FE fluctuations. Such fluctuations are indeed very large 
for the self-affine fractals that give rise to hot spots, and 
present a qualitative behavior similar to that of the maximum 
local FE in the vicinity of the hot spots. This is an interesting 
property that could be used in, e.g., PSTM studies to identify 
samples with the potential capability of yielding large optical 
excitations: even if no hot spots are found in the region being 
scanned, the calculated fluctuations of the resulting intensity 
map could provide a statistical account on the probability of 
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finding hot spots (simpler than calculating the total PDF for 
which much more data from a larger scanning area would be 
required). 

With regard to the quantitative aspects of SERS, the 
maximum local enhancement factors are still below those 
that could be deduced.-from experimental works on sin- 
gle molecule-jdetectiorffia and from approximate theoretical 
calculations OHj This, however, is not entirely surprising, 
given the differences in the type of SERS substrates being 
considered. In fact, the typical SERS spectroscopy enhance- 
ment factors are fairly similar to those found in this work. 

We would like to emphasize that our formulation is exact 
within the classical EM framework, at least as far as the lin- 
ear (direct) field enhancement factor is concerned. Thus, our 
calculations should provide a truthful picture of the linear op- 
tical response of self-affine fractal metal substrates. Further 
work is of course needed to test the perhaps naive assumption 
that the enhancement factor at the Raman-shifted frequency 
is identical to that obtained at the excitation frequency. Also, 
more work is required to study the effects of lower scaling 
cutoffs in the generation of the fractal surfaces, as these spa- 
tial frequencies might contribute to build up the enhancement 
factors. cj 

Finally, we mention that work involving rigorous calcula- 
tions of the kind presented here is also in progress for the 
study of the field enhancements produced on self-similar sub- 
strates, such as those found in colloidal aggregates,E3 and for 
the study surfaces covered by a monolayer of Raman active 
molecules (Langmuir-Blodgett films). 
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x (|j.m) 

FIG. 3. Surface EM field intensity for the region of the Ag fractal 
profile (also depicted) shown in Fig. ^| Solid curve: Electric field; 
dashed curve: Electric field, normal component; long-dashed curve: 
Electric field, tangential component; dot-dashed curve: Magnetic 
field. The inset zooms in the largest hot spot. 



(a) (b) 




FIG. 2. Near field intensity images (in a log 10 scale) resulting 
from the p-polarized scattering with 6q — 0°, A = 514.5 nm, 
and W = L/4cos#o, from a Ag fractal surface with D — 1.9, 
S = 514.5 nm, L = 10.29^m, and N p = 2000. The area shown is 
1 x 1/im 2 . (a) Electric field; (b) Electric field, z component (vertical); 
(c) Electric field, x component (horizontal); (d) Magnetic field. 



FIG. 4. Near electric field intensity image (in a log 10 scale) for 
the hot spot area shown in Fig. ^ for a Ag fractal surface with 
0o = 0°, D = 1.9, 8 = 514.5 nm, A = 514.5 nm, L = 10.29/xm, 
W = L/4cos6> , and N v = 2000. The area shown is 500x500 
nm 2 . 
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FIG. 5. Surface electric field intensity at the hot spot shown 
in Fig. ^ but for additional, slightly shifted incoming wavelengths 
A =476.9, 495.9, 539.1, and 563.6 nm. 



FIG. 6. Spectral dependence of the maximum local FE a max re- 
sulting from the p-polarized 
scattering with 9 = 0°, 5°, 10°, . . . , 50°, and W = L/4cos6> , 
from fractal metal surfaces with D = 1.9, consisting of N r = 60 
realizations of L = 10.29/im and jV p = 1600. (a) Electric field; (b) 
Electric field, normal component; (c) Electric field, tangential com- 
ponent; (d) Magnetic field. Circles, Ag; Squares, Au; Triangles, Cu. 
Filled symbols: 8 — 514.5 nm; Hollow symbols: 8 — 102.9 nm. 
Stars: water/ Ag, 8 = 514.5 nm. 
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FIG. 7. PDF of the p-polarized FE factor p(cr) for the 
electric (including both normal and tangential components) and 
magnetic fields resulting from the p-polarized scattering with 
0o = 0°, 5°, 10°, . . . , 50°, and W = L/4 cos O , from fractal metal 
surfaces consisting of iV r = 60 realizations of L = 10.29/im and 
iV p = 1600. Solid curve: D = 1.9 and 6 = 514.5 nm; dashed 
curve: D = 1.9 and 8 = 102.9 nm; and dotted curve: D — 1.2 and 
<5 = 102.9 nm. (a) A = 514.5 nm and (b) A = 1064 nm. 
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FIG. 8. Spectral dependence of the average surface FE (a) 
(excitation spectra) resulting from the p-polarized scattering with 
6 = 0° , 5° , 10° , . . . , 50° , and W = L/4 cos 6 , from fractal metal 
surfaces with D — 1.9, consisting of N r = 60 realizations of 
L = 10.29^m and N p = 1600. (a) Electric field; (b) Electric 
field, normal component; (c) Electric field, tangential component; 
(d) Magnetic field. Circles, Ag; Squares, Au; Triangles, Cu. Filled 
symbols: S = 514.5 nm; Hollow symbols: S — 102.9 nm. Stars: 
water/Ag, S = 514.5 nm. 
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FIG. 9. Same as in Fig. pmut for the absorption spectra at normal 
incidence. 
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